Brownian ratchet in a thermal bath driven by Coulomb friction 
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The rectification of unbiased fluctuations, also known as the ratchet effect, is normally obtained 
under statistical non-equilibrium conditions. Here we propose a new ratchet mechanism where a 
thermal bath solicits the random rotation of an asymmetric wheel, which is also subject to Coulomb 
friction due to solid-on-solid contacts. Numerical simulations and analytical calculations demon- 
strate a net drift induced by friction. If the thermal bath is replaced by a granular gas, the well 
known granular ratchet effect also intervenes, becoming dominant at high collision rates. For our 
chosen wheel shape the granular effect acts in the opposite direction with respect to the friction- 
induced torque, resulting in the inversion of the ratchet direction as the collision rate increases. We 
have realized a new granular ratchet experiment where both these ratchet effects are observed, as 
well as the predicted inversion at their crossover. Our discovery paves the way to the realization of 
micro and sub-micrometer Brownian motors in an equilibrium fluid, based purely upon nano- friction. 

PACS numbers: 02.50.Ey, 05.20.Dd, 81.05.Rm 
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From microscopic organisms to muscle fibres, from 
electric motors to power stations, the biosphere, our so- 
ciety and our lives critically depend on the conversion of 
energy to mechanical work. Thermodynamics provides 
precise and well established rules for energy conversion 
in macroscopic systems but these rules become blurred 
at small scales when thermal fluctuations play a decisive 
role Extracting work under such conditions requires 
subtle strategies radically different from those effective 
in the macroscopic world 0-0 • Within this framework, 
the theory of Brownian motors deals with the rectifica- 
tion of thermal fluctuations, a goal which can only be 
achieved in the presence of dissipation [|-13|. An in- 
teresting class of systems, where both dissipation and 
fluctuations are relevant, is represented by granular me- 
dia [l4, 15 1 . Indeed, interactions in a granular system are 
inherently dissipative, and because of its small number 
of constituents when compared with molecular gases or 
liquids, a granular fluid presents large fluctuations. The 
additional break of spatial symmetry is sufficient for a 
motor effect to be gen erated as demonstrated in a series 
of experiments [lj[3| and theoretical works 1^- 24 1 . 

In previous work the main source of dissipation was 
provided by the inelasticity of collisions, a property nor- 
mally not present at micro or nanometric scales. The 
remarkable result of our study is a new minimal model 
for a motor where energy is extracted from an equilibrium 
bath and dissipated only through Coulomb friction [2o| . 
Friction is therefore demonstrated to be an unexpect- 
edly efficient source of dissipation, that is able to rectify 



unbiased fluctuations also in the case of fully elastic col- 
lisions. Such a model can therefore be exploited in micro 
and nano apparatuses where friction is still present [26| . 
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FIG. 1: A) Sketch of the model, front view. B) Top view, 
with explanation of quantities used in the text for a generic- 
shaped rotator. C) Top view, with specific shapes used in the 
simulations and in the experiment. 



Our model, described pictorially in Fig. la, consists of: 
a wheel of mass M and inertia /, rotating with angular 
velocity u> around a fixed axis (say z). The wheel is 
immersed in an equilibrium fluid and collides with the 
molecules of mass to, and is subject to a viscous drag 
— Tyiscto and, most importantly, to a Coulomb friction 
torque —Ff r i C u on a(uj) (where a(x) is the sign function), 
due to solid-on-solid contacts within its support, e.g. a 
spherical bearing. The equation of motion for the angular 
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velocity u(t) of the wheel, therefore, reads 

u(t) = - j a u(t) - a[u(t)]A + i] coU (t) 



(1) 



where j a = T vlsc /I, A = Ff riction /I and rj coU (t) is the 
random force due to collisions with the molecules of the 
bath. The wheel is a cylinder parallel to the rotation axis 
z. Its base in the plane xy (shown in Fig. lb for a generic 
shape) can be symmetric or asymmetric for inversion of 
one of the two axis (x or y) . The two specific shapes taken 
in consideration here, one symmetric and other asymmet- 
ric, are drawn in Fig. lc. The velocities of the molecules 
are distributed according to the Maxwell-Boltzmann dis- 
tribution, with the only parameter being the "thermal" 
velocity vq = U (v 2 ) , where v is a component of the ve- 
locity vector on the xy plane. The molecular bath is 
also characterized by its number density n. The main 
parameter for the collision between the wheel and the 
molecules is the total cross section S. The collision rule, 
given in details in the Supplemental Material [l^, con- 
serves total angular momentum and may dissipate part 
of the total kinetic energy, according to the value of the 
restitution coefficient a G [0,1]. We will show that our 
model exhibits the ratchet effect even in the case of fully 
elastic collisions (a = 1), and even if the viscous force is 
removed — 0). The choice of a more general (pos- 
sibly inelastic) collision rule and the presence of a weak 
viscous damping is necessary to account for the results 
of the granular experiment described below. For consis- 
tency with the experiment, the viscous force (if present) 
is assumed to be small enough that 7 a |w| <C A for most 
of the values of ui. 

When the range of interactions with the molecules is 
short enough (as in the hard-core case), only two time- 
scales are relevant in the system: 1) the mean stop- 
ping time due to environmental dissipation, dominated 
by Coulomb friction, ta = ^"P*" ~ iFaj where 



pc 



denotes a post-collisional average, Ri = y/T/M is the 
radius of inertia and e = \JmjM ; 2) the mean free time 



between two collisions r c 
main control parameter 



We therefore use as 



(2) 

% /2vri? / A t c 

which is an estimate of the ratio of those two time- 
scales, as verified by simulations [13] ■ As already no- 
ticed [22j,|23|, when > 1 (t c < ta) the dynamics of 
the rotator is dominated by collisions and friction is neg- 
ligible (frequent collisions limit, denoted in the following 
by FCL); in the opposite limit fi~ x <C 1 (t c ^ ta, rare 
collisions limit, RCL) the rotator remains most of the 
time at rest and is rarely perturbed by collisions acting 
as independent random excitations. 

The complex behaviour of the model is simplified in 
the diluted limit, when Molecular Chaos can be as- 
sumed. With such an assumption, the probability den- 
sity function (pdf) p(cu,t) for the angular velocity is 



0.002 



0.001 



A 

G 
V 



-0.001 



-0.002 



#-# Right (a=0.4) 
A A Right (a-0.S 
Right (a=1) 
■■ h Symmetric (a-0.4) 
- - Theory 




1/P 



FIG. 2: Simulations under the assumption of Molecular 
Chaos. The rescaled average angular velocity (f2) is shown 
as a function of f3 . Theoretical expectations in the FCL 
and RCL are marked by dashed lines. The lower inset zooms 
in the RCL region. Simulations are performed using shapes 
and dimensions of Fig. lc, A = 10, 7 a = 0, Ri/e = 10 3 and 
Do — 100 in arbitrary units (varying nS to obtain a variation 
of/3" 1 ). 



fully described by the following linear Boltzmann equa- 
tion m m m 

d t p(ui, t) = d u [(Act(w) + 7 a w)p(cj, t)} + J\p, 0] (3a) 
J\p, 4>]= [ du?W{u\uf)p{uf, t) - p(w, t)/ c (w), (3b) 



where we introduce the rate W(u/|a;) for the transition 
uj — \ oj' and the velocity-dependent collision frequency 
/ c (w) = / duj'W(uj'\oj). The rate W(u)'\w), given ex- 
plictly for hard-core interactions in (2?]], depends on the 
velocity distribution of the gas particles, on the restitu- 
tion coefficient a, on the rotator cross section and on the 
density of the gas. 

A first insight into Eq. ([3]) is obtained by Direct Sim- 
ulation Monte Carlo [27], |29( , whose results are summa- 
rized in Fig. 2, always keeping j a = 0. The Figure shows 
the average velocity of the ratchet rescaled by the ideal 
"thermal" velocity, i.e. (Q) — J^-(cj), for several values of 
a and different shapes. Our main, unprecedented, result 
is the existence of an average drift, i.e. a motor effect, in 
the case of elastic collisions, provided that the shape is 
asymmetric (curve with diamond symbols). This effect 
is independent of the presence of the viscous damping, 
which very weakly affects the results of the simulation. 
In the elastic case, the average drift disappears for large 
where the effect of friction is washed out by highly 
frequent collisions and the system equilibrates with the 
bath. Remarkably, the ratchet effect starting from the 
RCL increases, in absolute value, when j3~ l increases, so 
that it must go through a maximum. We interpret such a 
maximum as a kind of stochastic resonance: the rotator 
switches from the "drift" state to the "rest" state on the 
time-scale ta and switches back to the "drift" state on 
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the time-scale r c . When the two scales synchronize with 
each other, the total time spent in the "drift" state is 
maximized, as well as its average velocity. 

As demonstrated in previous studies in the 

inelastic case (a < 1) a ratchet effect survives also in the 
FCL: interestingly, it takes different signs with respect 
to the RCL. Therefore, the crossover between these two 
limits requires the presence of an inversion point. 

It is possible to have an analytical account of the two 
opposite RCL and FCL limits [22|, |24j . When the mass 
of the rotator is large, e <C 1 the FCL is perturbatively 
reduced to a Brownian approximation and the average 
drift has already been computed [24] giving 
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the wheel Qsurf z 
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where the asymmetry of the rotator is represented by 
Afcl which is for a symmetric rotator; above we 
have used the shorthand notation for the uniform average 
along the perimeter (denoted as "surface" ) of the base of 
Isurf IT being the total perime- 
■jrp see Fig. lb for an explanation 
of symbols. This formula predicts zero net drift either 
with elastic collisions (a = 1) or with a symmetric rota- 
tor (Afcl — ), as expected from symmetry arguments. 
Most importantly, it predicts a constant value, as verified 
in the numerical simulations. This implies ~ vq for 
the dimensional angular velocity. 

The study of the RCL leads to remarkably different 
predictions. In such a limit, the dynamics after each 
collision event produces an increment of the angular po- 
sition of the rotator A6 which depends on the velocity v 
of the gas particle, precisely on its projection v = v • n, 

and on the point of impact represented by its curvilin- 

2 

ear abscissa s. The formula is A0(v,s) — ct^o)^ with 

2 

luq = — (1 + a) -gj i+Jg* ■ Following the calculations de- 
tailed in [27j, one obtains the formula for the rescaled 
average velocity of the ratchet 
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where Arcl = for symmetric shapes of the rotator. 
Eq. ([5]) shows that a non-zero drift is achieved for any 
value of the restitution coefficient: even in the (ideal) 
elastic case, Coulomb friction alone produces the desired 
ratchet effect provided that the shape of the rotator is not 
symmetric, i.e. that Arcl 7^ 0. Note that the limit of 
vanishing dry friction (A — > 0) is singular in formula (|5a[) , 
since in the absence of dissipation between collisions the 
stopping time becomes infinite, ta — ► 00, and the as- 
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FIG. 3: Real experiment. A) Experimental setup. B) The 
angular position of the rotator as a function of time, high- 
lighting the ratchet effect for the asymmetric rotator. The 
inset shows an enlargement (ten minutes) extracted from a 
asymmetric (right) experiment. For these experiments the 
value of maximum acceleration, normalized by gravity, is 13. 



sumption of "rare collisions" breaks down. Equally inter- 
esting, the shape factor Arcl, determining the intensity 
and drift direction in the RCL, can take opposite sign 
with respect to the shape factor Afcl in the FCL for- 
mula. This is precisely the case for our chosen shape, 
see Fig. lc. Moreover, the magnitude of the drift is pre- 
dicted to increase with /3 _1 as seen in the numerical sim- 
ulations for small /3 _1 . This corresponds to ~ Vq. 
Both the predictions for the RCL and for the FCL are 
superimposed on the results (elastic and inelastic) of the 
numerical simulations in Fig. 2, demonstrating excellent 
agreement in their respective limits. We remark that if 
friction is removed (A = 0) the only ratchet effect is the 
one predicted in Eq. (j4|) , i.e. no inversion is observed. 

In order to obtain the first experimental evidence of 
this newly discovered ratchet effect, we have built a 
macroscopic realization of our model, i.e. a setup where 
the thermal bath is replaced by a fluidized granular gas. 
In such a setup the collisions are unavoidably inelastic: 
nonetheless, by tuning the collision frequency, it is possi- 
ble to disentangle the two ratchet mechanisms which act 
in opposite directions, and so demonstrate the newly dis- 
covered effect induced by Coulomb friction. Our setup 
consists of a rotator vertically suspended in a granular 
medium (see Fig. 3a) maintained by an electrodynamic 
shaker in a (roughly homogeneous) stationary gaseous 



regime 
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The statistics of the velocities of the 



grains, on the rotation plane, has been verified to be 
indistinguishable from a Maxwell-Boltzmann distribu- 
tion [271 ] . The shaker performs a vertical sinusoidal os- 
cillation at a frequency of 53 Hz, while the amplitude 
is varied to explore different regimes of the system. We 
stress that the rotator is not in direct contact with the 
shaker, it only collides with the flying grains. Its motion 
is recorded by an angular encoder which also supports it 
through two precision spherical bearings. Two rotators 
have been realized to reproduce the two different shapes 
of Fig. lc. The asymmetry of the latter can be inverted 
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FIG. 4: Real experiment. Plot of the rescaled angular velocity 
of the rotator, (D.) = -^(uj), averaged on an 8-hours run, as 
a function of /3 _1 which estimates the ratio of time-scales 
The main plot shows the region where a maximum and 
a current inversion are observed for the asymmetric rotator. 
The inset sums all the results up, highlighting the behaviour 
for large /3 _1 . Experiments are performed with maximum 
shaker acceleration varying in the range from 5 to 20 in gravity 
units. 



(from left to right hand, and vice versa) simply by turn- 
ing the rotator upside-down. The amplitude of the shaker 
oscillation is varied to span a range of maximum accel- 
eration (in units of gravity acceleration) between 5 to 
20, resulting in a granular thermal velocity vq ranging 
between ~ 100 mm/s and ~ 600 mm/s. Details of the 
experimental setup and measurement of the parameters 
are in [13]. 

In Fig. 3b, we provide some runs demonstrating the 
ratchet effect for the asymmetric rotator. A net drift 
is evident when a chiral rotator is used. Turning the 
rotator upside down results in a reversed direction of ro- 
tation. The symmetric rotator makes a diffusive motion 
with only a small drift revealing the presence of some 
bias due to imperfections in the rotational symmetry of 
the setup. Such a bias affects also the curves pertaining 
to the asymmetric rotator and it does not appear to de- 
pend on the rotator's direction of asymmetry, or absence 
thereof. The typical instantaneous velocity of rotators 
goes from ~ 10 _1 to ~ 1 rad/s, much larger than the 
typical ratchet drift, which goes from - 10~ 4 to 10" 2 
rad/s. Such a low signal-to- noise ratio makes it impossi- 
ble to determine the true behaviour of the rotator if one 
monitors its position only for short times (e.g. less than 
a few hours), as evidenced in the inset of Fig. 3b. 

In Fig. 4, we show the average adimensional angular 
velocity (f2) = ^j- (oj) for a set of experiments with asym- 
metric and symmetric rotators. The behaviour with fi 
strongly resembles that observed in our numerical simula- 
tions (Fig. 2). At /3 _1 < 1, we measure |(fi)| increasing 
with P~ l , followed by a maximum in the proximity of 



~ 1 and then by an inversion of direction of motion. 
At large , | (fl) | increases and finally reaches a plateau 
(see the inset of Fig. 4) . Even with some quantitative dis- 
crepancies, we can claim that our experiment reproduces 
very well the qualitative phenomenology of the model, in- 
cluding the resonant maximum and the inversion point, 
which are both evidence of the presence of the friction- 
induced ratchet mechanism. We believe that the quanti- 
tative differences (the real ratchet is faster roughly by a 
factor 2 in the RCL and a factor 5 in the FCL) can be 
imputed to the many assumptions present in the model, 
the most important being molecular chaos and spatial ho- 
mogeneity, hardly controlled in the experiment: indeed 
non-equilibrium correlations may well become important 
at high collision frequencies [3l| . 

To conclude, our main discovery is the existence of 
a minimal ratchet model made of two simple ingredi- 
ents: a wheel subject to Coulomb friction and a bath 
at thermodynamic equilibrium. Such a model appears 
even simpler than the classical Feynmann-Smoluchowsky 
model [9j . The observation in the laboratory of a maxi- 
mum and an inversion of the ratchet velocity (Fig. 4) , due 
to the crossover from the inelasticity-dominated (FCL) 
to the friction-dominated (RCL) regime, is a strong ex- 
perimental demonstration of the efficiency of this effect. 
We wish to remark that, in all previous experimental 
and theoretical work on friction-driven ratchets (32T - [3o| . 
the energy injection was provided by mechanisms differ- 
ent from an equilibrium bath: our proposal is the first 
which can be realized at the micro- and nano-scale in an 
equilibrium fluid, that is, without the application of any 
external field. 

We would like to thank A. Vulpiani for useful com- 
ments and MD. Deen for technical support. The au- 
thors acknowledge the support of the Italian MIUR un- 
der the grants: FIRB-IDEAS n. RBID08Z9JE, FIRBs 
n. RBFR081IUK and n. RBFR08M3P4, and PRIN n. 
2009PYYZM5. 



SUPPLEMENTAL MATERIAL 
Details of the theory 

We recall here some of the main quantities of our model 
(refer to Figure 1 for a visual explanation of symbols). A 
rigid body of momentum of inertia / and total mass M 
is bound to rotate around a fixed axis (say z) and is sus- 
pended in a molecular or granular fluid. The body is 
constituted by the set of material points with cartesian 
coordinates {x, y, z} with z g [0, h] (where h is the height 
of the cylinder) and \/ x 1 + y 2 < r(s) for each s <E [0, S] 
where s is the curvilinear abscissa, r(s) is the curve de- 
limiting a section of the solid in the xy plane, and S is 
the perimeter of the section. The fluid surrounding the 
rotator has number density n and is made of identical 
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spheres of mass m. We denote by us the angular velocity 
of the rotator, by 9 its angular position, and by v the 
velocity of a molecule of the fluid. We also denote by 
p = nh the two-dimensional projection of density, which 
is the only one which matters in our problem. Note that 
pS = nYj if £ (as in the Letter) is the total scattering 
cross section of the rotator against the molecules of the 
bath. 

Inelastic collisions The fluid interact with the rotator 
by means of inelastic collisions, which change the velocity 
of the rotator from u> (before the interaction) to u/ (after 
the interaction) and that of the colliding particle from v 
to v'. In the case of hard-core collisions the scattering 
event follows the rules 



u = oj + (1 + a) 



(V — v) • h ge 2 



Rj 



1 + eV 



, (V - v) • n „ 
y > = v + 1 + a y ' n 
1 + e z g z 



(6a) 
(6b) 



where a € [0, 1] is the restitution coefficient (a = 1 cor- 
responds to elastic collisions), V = cuz x r is the linear 
velocity of the rotator at the point of impact r, n is the 
unit vector perpendicular to the surface at that point, 
and finally g = with t — z x h which is the unit 
vector tangent to the surface at the point of impact. 

Equations ([5]) guarantee that total angular momentum 
L'z = mr x v+iwz is conserved, that relative velocity pro- 
jected on the collision unit vector is reflected and rescaled 
by the restitution coefficient, (V — v')-n = — a(V— v)-n, 
and finally that the kinetic energy K — ?r|v| 2 
changes as 



1 ' , , 2 



K'-K = - {l - a2 \(y- 
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(7) 



A few relations in cartesian coordinates may be useful: 
V = (— wr y , u)r x ) and t — (—n y ,n x )- It is also useful to 
realize that V • h = —to Rig. 

Transition rates for a Gaussian gas The transition rate 
for the collisional Markov process reads [24j 



W(u'\u)= P S J^-J dv0(v)0[(V-v) 
|(V-v) ■ h\5[u>' - w - Acj(s)], 



n x 



Aui(s) = (1 + a) 



[V( S ) 



- U) 

vl • n 



Ri{s) l + e 2 5 (s) 2 ' 



(8a) 
(8b) 
(8c) 



where 4>(y) is the pdf for the gas particle velocities and 
the Heaviside step function 0[(V — v) • n] enforces the 
kinematic condition necessary for impact. 

Assuming that the velocities of the molecules of the 
surrounding fluid obey the Maxwell-Boltzmann statistics 
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one gets an explicit expression for the 



transition rates 
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Average velocity in the RCL When /3 _1 <C 1 (ta <C t c , 
denoted as Rare Collisions Limit, RCL), the dynamics is 
dominated by the rotator at zero velocity with random 
and independent perturbations due to sparse collisions 
with gas particles [22[. In this case, the dynamics after 
each collision event produces an increment of the angular 
position of the rotator A8 which depends on the velocity 
v of the gas particle, precisely on its projection v = v • n, 
and on the point of impact represented for instance by 
its curvilinear abscissa s. The formula in the limit of 
negligible viscous damping (j a —> 0) is 



uj = -(1 + a) 



v e-g 



Ri 1 + e 2 g 2 ' 



(10) 
(11) 



As a consequence, the average velocity of the ratchet 
reads [22 1 



/ x PS f i i i , / \ n > s f ds , s v 2 (1 + a) 2 e 4 g 2 

<"> = £a / dv\v\<Kv)Q(-v) I -cT(g)— \^ J 2g2 * 2 

(12) 
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(l + eV)V W 



where we have used the shorthand notation for the aver- 
age along the perimeter of the rotator's shape () S urf — 
J 4p If one considers the expression for the control pa- 
rameter /3 _1 we can finally write for the adimensional 
angular velocity f2 = — u 



(n) = V^(l + a) 2 (3- 1 e 2 (^j i 



v(g)g 2 



+ ( 2 9 2 ?/ surf 



(13) 



Details of the simulations 

We have simulated the model described in equa- 
tions (2) of the Letter, through a suitable adaptation of 
the Direct Simulation Monte Carlo method (DSMC) [29| . 
The DSMC is devised to solve numerically a Boltzmann 
equation, therefore enforcing the Molecular Chaos as- 
sumption. In our specific problem, however, the proce- 
dure is drastically simplified, since only one particle (the 
rotator) is represented in the simulation, through its an- 
gular velocity and position, uj and 9 respectively. The 
surrounding gas is represented by its constant velocity 
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distribution <fc(v) (assumed Gaussian with variance Vq), 
unaffected by collisions. The dynamics of these variables 
advances by a series of constant and small time steps of 
length St taken smaller than all characteristic time-scales 
in the problem. At every time step a free streaming up- 
date and a collisional update are performed. The free 
streaming corresponds to the evolution of u> and 6 from 
t to t + St in the absence of any interrupting collisions. 
In the collisional update, a collision is performed with 
the probability dictated by the Boltzmann equation (Eq. 
(3) of the Letter). The correct probability is sampled 
through a Monte Carlo procedure, where a tentative col- 
lision is proposed, by choosing the point of impact (at 
random and uniformly) along the surface of the rotator 
and by extracting the velocity v with probability </>(v). 
The tentative collision is realized with a probability equal 
to p c = (V — v) • hpSSt (note that St is chosen much 
smaller than r c and this guarantees that p c <C 1), lead- 
ing to an update of lo by the collision rule. Otherwise no 
collisions occur. Note that the choice of the point of im- 
pact (which determines also n) can reproduce the precise 
shape of the experimental rotator (symmetric or chiral). 

Control parameter /3 _1 as a good estimate of the ratio 
of typical times. In Figure 5 we report the ratio between 
the collision times r c and frictional stop times ta, ob- 
served in the simulations in different regimes, in order to 
assess the fairness of estimate (5~ x « ta/t c . 
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FIG. 5: Assessment of the estimate for ta/t c . The ratio 
between ta/t c and obtained in numerical simulations is 
shown in the inset. The estimate is useful in experiments, 
where measurements of ta and r c are not reliable. It is seen 
that the estimate is a good approximation of the time-scale 
ratio where it is close to 1, while it slightly deviates from it in 
the RCL and FCL, by a factor ~ 1.25 and ~ 0.7 respectively. 



Experimental set-up. 

We assume that the evolution of the angular velocity 
lu of the rotator obeys an equation of the form 



Iu>(t) = -Ff riction a[uj(t)] 



c ij(t)+F coU (t), (14) 



where / is the momentum of inertia, Ff r i Ct i on is the dry 
friction torque, a the sign function and T V i SC the air vis- 
cous drag. The term F co u — Ir\ co ii contains the driving 
torque, that is the angular momentum randomly trans- 
ferred from the beads to the rotator, and is analysed in 
detail in Section 2. We refer to Figure 3a of the paper for 
the main features of our set-up. Dissipation coefficients 
rescaled by inertia are 7 a = T msc /I and A = Ff riction /I. 

Technical details about the experimental setup. The 
granular medium, made of 50 polyoxymethylene (POM) 
spheres (radius r — 3 mm and mass m = 0.15 g), 
is housed in a polymethyl-methacrylate (PMMA) cylin- 
der (diameter 90 mm) with a conical-shaped floor. A 
removable cap encloses a miniaturized angular encoder 
(model AEDA-3300 by Avago Technologies). The en- 
coder, which also supports the rotator, provides high res- 
olution measurements (up to 80,000 division/revolution 
at the maximum rate of 20 kHz) of the rotator posi- 
tion. The system is vibrated by an electrodynamic shaker 
(model V450 by LDS Test & Measurement) fed by a sinu- 
soidal excitation. An accelerometer measures the actual 
acceleration induced to the system. A high-speed camera 
(EoSens CL by Mikrotron) tracks single beads in order 
to measure their velocity. The two rotators are made of 
PMMA, have height h = 15 mm and are distinguished 
by their different section on the rotation plane. The mo- 
mentum of inertia of the rotator / comprises different 
parts. The probe is attached to the angular encoder by 
means of a steel rotation axis 50 mm high and 3 mm 
thick. The mass of the whole rotator (axis and probe) is 
M = 5.21 g for the asymmetric probe and M = 6.49 g 
for the symmetric one. The total momentum of iner- 
tia is / = 135 gmm 2 and / = 353 gmm 2 for the two 
types, respectively. In both cases the inertia of the axis 
is a few hundredths of the total one. The shape factors 
of the asymmetric rotator, relevant for its ratchet effect, 
T -te£-\ - 0.0013 and ^ 



the RCL and FCL limits respectively (see text for defini- 
tions). Both the shape factors vanish for the symmetric 
rotator. 

Dry and viscous friction. The only source of dry fric- 
tion comes from the two ball bearings inside the encoder 
while viscous friction is due to the air surrounding the ro- 
tator. We measure the dry and viscous friction during a 
standard experimental run (eight hours) , but considering 
only periods where the rotator is subject to these forces 
and no other: lu — ~-f a uj^fA (where the upper sign holds 
for lj > and the lower for u> < 0), i.e. excluding the time 
immediately about collisions and periods of rest. Angular 
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velocity and acceleration are obtained from the angular 
position (in function of time) by differentiation. A bidi- 
mensional histogram is calculated and a peak-finding pro- 
cedure has been applied in order to obtain the most prob- 
able trajectories in (lo,u>) space. Finally, the most prob- 
able trajectory is fitted with the above linear equations, 
yielding the values Ff ricUon — (9.9 ±0.7) x f0 3 gmm 2 /s 2 
and T visc = (1.6 ± 0.2) x 10 3 gmm 2 /s. It is inter- 
esting to notice that viscous friction becomes impor- 
tant for velocities larger than the threshold velocity 
uth = Ff„ ction /T vlsc w 6.2 s~\ separating the friction- 
dominated regime u <C uith horn the viscosity-dominated 
one oj ^> LOth- 

Restitution coefficients. The beads-rotator restitution 
coefficient a has been measured launching a single bead 
against the rotator while recording the rotator position at 
high sampling rate (1 kHz). For these measurements the 
top of the apparatus is removed and put with the rota- 
tion axis parallel to the floor. A bead falls from height h 
and hits the rotator with velocity v — yj2a g h (here a g in- 
dicates the gravity acceleration). The high-speed camera 
has been used to monitor the exact distance x of the im- 
pact point from the rotation axis. We calculate the resti- 
tution coefficient using the collision rule (see Section 2 
below) adapted for this particular configuration in which 
the rotator is at rest before the collision (51 = V = 0) 
and g ~ xjRj. In this circumstances the velocity of the 
rotator after a collision is a/ = (l + a)vx/(I/m+x 2 ). Re- 
peated measurements on the symmetric rotator gave the 
following results: a — 0.83 ± 0.16 (we recall that a = 1 
for elastic collisions). 

Granular gas velocity statistics. We discuss here the 
methods used to characterize the velocity statistics of the 
granular gas. A fast camera (EoSens CL by Mikrotron) 
is placed above the PMMA container to catch horizontal 
components of motion (y x ,v y ) of the polyoxymethylene 
(POM) spheres constituting the granular gas. Focus of 
lens is adjusted to the plane at half height of the probe. 
Lighting is provided by a led-dome which diffuses light 
on the system reducing shadows and reflections. Particle 
tracking is enhanced by marking a few tracers: 3 POM 
spheres are white, while all the other particles are black. 
Pictures are taken at 250 frames per second: this is veri- 
fied to be the optimal compromise between too large de- 
lays, which prevent tracking of ballistic trajectories, and 
too small ones, which produce "false movement" induced 
by noise due to finite sensitivity of the camera resolution 
and the centring algorithm. Uncertainty in the determi- 
nation of the centre of mass of tracer spheres is estimated 
to be ~ 0.05 mm. The details of velocity measurements 
and error estimates are reported in [36j |. In Figure 6 we 
display a typical snapshot of the system from the camera 
(see also the accompanying Movie). 

For each choice of shaking parameters we have com- 
puted the velocity histogram from a series of 2 x 10 4 
contiguous frames. Examples of those histograms are re- 



FIG. 6: A snapshot of the system taken from the fast camera 
(shaker is vibrating at 53 Hz with rescaled maximum acceler- 
ation r = 8). 

ported in Figure 7, left frame. A Gaussian fit is seen 
to be a good approximation and provides the value of 
v = v (v 2 ) where v is v x or v y (isotropy is always veri- 
fied) . In the right frame of Figure 7 we report the values 
of vq as function of the rescaled maximum acceleration 
r = a '^° x where a g is gravity acceleration and the z s (t) 
position of the shaker's vibrating head follows the law 
z s (t) = &™j?2 sin(27r/t) with / being the vibration fre- 
quency. 

Rotator velocity signal. An example of the velocity 
signal u)(t) in two different regimes (frequent and rare 
collisions) is shown in Figure 8, left frame, for the case 
of the symmetric rotator. In the right frame of the 
same Figure, we show the angular velocity autocorrela- 
tion C(t) = (uj(t)uj(Q)) — (uj) 2 versus time. Note that the 
signal average is much smaller than its standard devia- 
tions and, hence, it is negligible in C(t). 

Here, it is interesting to comment about the smaller 
characteristic time associated to the system with rare col- 
lisions (red curve) , with respect to that with frequent col- 
lisions (black curve) : this is a consequence of the smaller 
typical velocity in the rare collisions configuration, which 
induces faster relaxations through Coulomb friction dis- 
sipation (we recall that ta ~ Similar results are 
obtained with the asymmetric rotator. 
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FIG. 7: Left: Histograms of velocities (v x ) of iV = 50 granu- 
lar gas particles for different values of the rescaled maximum 
acceleration V at vibration frequency 53 Hz. Gaussian fits are 
superimposed. Right: the values of vo (standard deviation) 
versus V. 
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FIG. 8: Left: Examples of symmetrical rotator's angular 
velocity signal u)(t) as acquired from the angular encoder at 
200 Hz in two different regimes (frequent collisions, top, and 
rare collisions, bottom). Right: rescaled angular velocity au- 
tocorrelation C(t)/C(0) versus time t for the two cases shown 
in the left frame. 
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